Sharp error estimates for a discretisation of the 1D convection/diffusion equation with Dirac initial data

This paper derives sharp l$\infty$ and l1 estimates of the error arising from an explicit approximation of the constant coefficient 1D convection/diffusion equation with Dirac initial data. The analysis embeds the discrete equations within a semi-discrete system of equations which can be solved by Fourier analysis. The error estimates are then obtained through asymptotic approximation of the integrals resulting from the inverse Fourier transform. this research is motivated by the desire to prove convergence of approximations to adjoint partial differential equations

Aerodynamic design optimisation for complex geometries using unstructured grids

These lecture notes, prepared for the 1997 VKI Lecture Course on Inverse Design, discuss the use of unstructured grid CFD methods in the design of complex aeronautical geometries. The emphasis is on gradient-based optimisation approaches. The evaluation of approximate and exact linear sensitivities is described, as are different ways of formulating the adjoint equations to greatly reduce the computational cost when dealing with large numbers of design parameters. \ud \ud The current state-of-the-art is illustrated by two examples from turbomachinery and aircraft design

Monte Carlo evaluation of sensitivities in computational finance

In computational finance, Monte Carlo simulation is used to compute the correct prices for financial options. More important, however, is the ability to compute the so-called "Greeks'', the first and second order derivatives of the prices with respect to input parameters such as the current asset price, interest rate and level of volatility.\ud \ud This paper discusses the three main approaches to computing Greeks: finite difference, likelihood ratio method (LRM) and pathwise sensitivity calculation. The last of these has an adjoint implementation with a computational cost which is independent of the number of first derivatives to be calculated. We explain how the practical development of adjoint codes is greatly assisted by using Algorithmic Differentiation, and in particular discuss the performance achieved by the FADBAD++ software package which is based on templates and operator overloading within C++.\ud \ud The pathwise approach is not applicable when the financial payoff function is not differentiable, and even when the payoff is differentiable, the use of scripting in real-world implementations means it can be very difficult in practice to evaluate the derivative of very complex financial products. A new idea is presented to address these limitations by combining the adjoint pathwise approach for the stochastic path evolution with LRM for the payoff evaluation

An extended collection of matrix derivative results for forward and reverse mode automatic differentiation

This paper collects together a number of matrix derivative results which are very useful in forward and reverse mode algorithmic differentiation (AD). It highlights in particular the remarkable contribution of a 1948 paper by Dwyer and Macphail which derives the linear and adjoint sensitivities of a matrix product, inverse and determinant, and a number of related results motivated by applications in multivariate analysis in statistics.\ud \ud This is an extended version of a paper which will appear in the proceedings of AD2008, the 5th International Conference on Automatic Differentiation

Adjoint recovery of superconvergent functionals from PDE approximations

Motivated by applications in computational fluid dynamics, a method is presented for obtaining estimates of integral functionals, such as lift or drag, that have twice the order of accuracy of the computed flow solution on which they are based. This is achieved through error analysis that uses an adjoint PDE to relate the local errors in approximating the flow solution to the corresponding global errors in the functional of interest. Numerical evaluation of the local residual error together with an approximate solution to the adjoint equations may thus be combined to produce a correction for the computed functional value that yields the desired improvement in accuracy. Numerical results are presented for the Poisson equation in one and two dimensions and for the nonlinear quasi-one-dimensional Euler equations. The theory is equally applicable to nonlinear equations in complex multi-dimensional domains and holds great promise for use in a range of engineering disciplines in which a few integral quantities are a key output of numerical approximations

Convergence analysis of Crank-Nicolson and Rannacher time-marching

This paper presents a convergence analysis of Crank-Nicolson and Rannacher time-marching methods which are often used in finite difference discretisations of the Black-Scholes equations. Particular attention is paid to the important role of Rannacher's startup procedure, in which one or more initial timesteps use Backward Euler timestepping, to achieve second order convergence for approximations of the first and second derivatives. Numerical results confirm the sharpness of the error analysis which is based on asymptotic analysis of the behaviour of the Fourier transform. The relevance to Black-Scholes applications is discussed in detail, with numerical results supporting recommendations on how to maximise the accuracy for a given computational cost

Computing Greeks using multilevel path simulation

We investigate the extension of the multilevel Monte Carlo method [2, 3] to the calculation of Greeks. The pathwise sensitivity analysis [5] differentiates the path evolution and effectively reduces the smoothness of the payoff. This leads to new challenges: the use of naive algorithms is often impossible because of the inapplicability of pathwise sensitivities to discontinuous payoffs.\ud \ud These challenges can be addressed in three different ways: payoff smoothing using conditional expectations of the payoff before maturity [5]; an approximation of the above technique using path splitting for the final timestep [1]; the use of a hybrid combination of pathwise sensitivity and the Likelihood Ratio Method [4]. We discuss the strengths and weaknesses of these alternatives in different multilevel Monte Carlo settings

Sharp error estimates for discretisations of the 1D convection/diffusion equation with Dirac initial data

This paper derives sharp estimates of the error arising from explicit and implicit approximations of the constant coefficient 1D convection/diffusion equation with Dirac initial data. The error analysis is based on Fourier analysis and asymptotic approximation of the integrals resulting from the inverse Fourier transform. This research is motivated by applications in computational finance and the desire to prove convergence of approximations to adjoint partial differential equations

Smoking Adjoints: fast evaluation of Greeks in Monte Carlo calculations

This paper presents an adjoint method to accelerate the calculation of Greeks by Monte Carlo simulation. The method calculates price sensitivities along each path; but in contrast to a forward pathwise calculation, it works backward recursively using adjoint variables. Along each path, the forward and adjoint implementations produce the same values, but the adjoint method rearranges the calculations to generate potential computational savings. The adjoint method outperforms a forward implementation in calculating the sensitivities of a small number of outputs to a large number of inputs. This applies, for example, in estimating the sensitivities of an interest rate derivatives book to multiple points along an initial forward curve or the sensitivities of an equity derivatives book to multiple points on a volatility surface. We illustrate the application of the method in the setting of the LIBOR market model. Numerical results confirm that the computational advantage of the adjoint method grows in proportion to the number of initial forward rates